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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  {AEDC),  Air  Force  Systems  Command  (AFSC).  The  results  of  the  research  were 
obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup  Corporation  Company),  operating 
contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project 
Number  V32I-P2A.  The  Air  Force  project  manager  was  Captain  Stanislaus  L.  Ludwig  (CF), 
AEDC/DOT.  The  data  analysis  was  completed  on  August  18,  1978,  and  the  manuscript  was 
submitted  for  publication  on  March  30,  1979. 
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1.0  INTRODUCTION 

There  is  a  continuing  need  for  measurement  capabilities  of  particulates  encountered  in 
ground-based  simulation  facilities.  These  particles  may  be  generated  by  the  source  of  the 
simulation  flow  field  or,  in  some  cases,  by  the  object,  such  as  a  rocket  engine,  which  is  being 
tested.  The  characterization  provided  by  the  diagnostics  can  be  quite  varied,  ranging  from 
mere  determination  of  the  existence  of  the  particulates  to  the  more  complete  specification  of 
particle  size  distribution  function  and  identification  of  particulate  species  and  the  three- 
dimensional  velocity  distribution  function.  For  ground-based  simulation  facilities  as  diverse 
as  those  of  the  Arnold  Engineering  Development  Center  (AEDC),  the  range  of  particle  sizes 
and  densities,  particulate  species,  flow  velocities,  required  time  and  space  resolution,  and 
spatial  extent  of  the  particulates  obviously  provides  multidimensional  diagnostics 
requirements.  Therefore,  it  should  not  be  surprising  that  no  one  measurement  approach  is 
adequate  for  the  spectrum  of  testing  requirements.  Simple  examples  of  the  restriction  of 
application  regimes  of  various  techniques  are  a)  the  limitation  of  current  holographic 
capability  by  particles  of  super-micron  sizes;  and  b)  the  limitation  to  low  particulate 
densities  of  techniques  which  require  the  presence  of  no  more  than  one  to  two  particles 
within  the  focal  volume  of  an  injected,  focused  laser  beam.  It  cannot  be  too  strongly 
emphasized  that  the  regimes  of  application  of  many  techniques  become  oppressively 
constraining  when  the  desire  for  short  time  resolution  in  pulsed  flow  facilities  is  introduced. 
Perspective  of  pulsed  flow  field-particulate  diagnostics  requirements  is  afforded  by  noting 
that  a  flow-field  duration  of  100  msec  and  a  desired  time  resolution  on  the  order  of  1  msec 
are  not  uncommon  test  parameters  either  for  ground-based  test  facilities  or  for  pulsed  liquid 
rocket  engine  plume  exhausts.  If  one  includes  certain  solid  rocket  engines  in  this 
consideration,  the  requirements  increase  in  stringency  by  an  order  of  magnitude  (i.e.,  to 
10-msec  duration  and  100-jtsec  resolution).  Clearly,  nonroutine  diagnostic  techniques  are 
required. 

Presented  in  this  work  is  the  description  of  a  particulate  visualization  technique  which 
was  employed  recently  for  characterization  of  a  nominal  100-msec-duration  flow  field  and 
for  which  time  resolution  on  the  order  of  1  msec  was  achieved.  As  will  be  presented  in  detail, 
the  application  of  interest  was  the  visualization  of  the  spatial  extent  of  particulates  injected 
into  a  coflowing  hypersonic  flow  field  which  was  provided  by  the  AEDC  von  Karman  Gas 
Dynamics  Facility  (VKF)  arc-heated  Hypervelocity  Wind  Tunnel  (F).  The  spatial  extent  of 
the  injected  particulates  w<as  relatively  large,  and  qualitative  knowledge  of  the  density 
distribution  was  desired.  Information  on  the  magnitude  of  heating  and  the  subsequent 
radiation  of  the  injected  particles  following  traversal  of  a  shock  system  was  also  desired. 

The  approach  taken  to  satisfy  these  test  requirements  was  the  injection  into  the  flow  field 
of  a  laser  beam  which  was  expanded  into  a  wide  plane  sheet  configuration.  The  interaction 
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of  this  beam  with  the  particulates  produced  Mie  scattered  radiation  by  which  the  scene  was 
viewed  and  which  was  detected  with  an  image  intensifier.  To  achieve  the  desired  time 
resolution  required  that  the  output  of  the  intensifier  be  viewed  with  a  rapid-framing  camera. 
Modulation  of  the  laser  beam  at  the  100-percent  level  provided  images  produced  by  the 
particulate  emission  alone,  and  these  data  yielded  information  regarding  particulate  heating. 
Data  reduction  and  analysis  were  achieved  using  computer  analysis  of  the  fast-scan 
photographs.  A  more  complete  description  of  the  system,  the  application  of  interest,  and  the 
results  is  presented  in  the  following  section.  Finally,  an  order-of-magnitude  comparison  of 
the  sensitivity  of  the  image  intensifier  detection  system  to  the  sensitivity  of  the  more 
conventional  photographic  method  is  provided. 

2.0  APPLICATION  CONFIGURATION  AND  THEORY 

2.1  CONFIGURATION  OF  APPLICATION 

The  particulate-laden  flow  field  which  is  to  be  considered  is  partially  sketched  in  Fig.  1, 
which  shows  the  source  from  which  the  particles  are  injected  into  a  coflowing  hypersonic 
flowfield.  (The  gas  dynamic  shock  systems  are  not  shown  in  Fig.  1.)  The  laser  beam  is 
expanded  to  illuminate  a  substantial  axial  length  of  the  particle  plume  exhaust.  The  Mie 
scattered  radiation  is  collected  and  detected  by  the  image  intensifier  optical  system. 

Prior  to  the  application  of  this  system  to  the  visualization  task,  as  shown  in  Fig.  1,  a 
feasibility  calculation  was  performed  with  the  goal  of  determining  the  photon  rate  of 
irradiation  of  a  picture  element,  or  pixel,  of  the  image  intensifier  tube.  The  results  of  this 
calculation  can  be  compared  with  the  calculated  performance  of  a  system  which  employs 
photographic  film  detection  rather  than  intensification  followed  by  photography.  The 
calculation  presented  in  the  following  sections  proceeds  thus:  a)  description  of  single-particle 
Mie  scattering;  b)  generalization  of  the  single-particle  theory  to  the  situation  as  shown  in 
Fig.  1;  and  c)  calculation  of  the  response  of  a  photographic  system. 

2.2  MIE  SCATTERING  THEORY  FOR  LOCALIZED  PARTICLES 

The  incoherent  scattering  of  an  incident  plane  wave  by  a  single  dielectric  sphere  is 
described  exactly  by  the  Mie  scattering  theory,  which  was  developed  in  the  first  decade  of  the 
20th  century  through  the  collective  efforts  of  Lorentz,  Mie,  Debye,  and  others.  A 
comprehensive  derivation  of  Mie  scattering  theory  is  given  in  Ref.  1.  Practical  application  of 
Mie  scattering  to  problems  of  current  interest  has  been  made  possible  comparatively  recently 
by  the  employment  of  modern  digital  computing  techniques,  which  are  exemplified  by  Ref. 
2.  For  the  purposes  of  the  calculations  presented  in  this  report  the  authors  use  a  computer 
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program  described  in  Ref.  3  and  refer  the  reader  to  the  literature  cited  therein  for  details  of 
the  calculation  methodology. 

Reference  to  Fig.  2  will  aid  the  reader  in  understanding  the  general  attributes  of  Mie 
scattering  from  a  single  particle.  A  plane  wave  of  wavelength  \  and  an  electric  field  polarized 
along  the  x  axis  are  incident  on  a  particle  of  complex  refractive  index  17  =  rjr  -  ix,  which  is 
located  at  the  origin.  The  direction  of  propagation  of  the  incident  wave  is  along  the  positive 
z  axis;  i.e.,  the  particle  is  illuminated  from  below.  Upon  interaction  with  the  particle  the 
light  is  scattered,  and  if  x  *  0,  a  portion  of  its  energy  is  lost  to  absorption  by  the  particle. 
For  an  observer  whose  coordinates  are  (r,  6,  0),  the  intensity  of  the  light  scattered  without 
change  in  wavelength  is  given  by 

lo 

I  {r,  6,  o,  rj,  xl  -  -  jjj  <x.  0,  if)  sin20  +  i9  (x.  0  tj )  cos2<iJ  (11  (1) 


w'here 

Io  is  the  incident  light  intensity  lenetgv  flux  pci  unit  time  per  umL  circa), 
k  -  2tt'X  is  the  wave  number  of  incident  light, 
r  is  the  (hf-tuncc  from  the  particle  to  the  observer. 

\  -  rrD-,\  is  the  particle  size  paiamclct  based  oil  dimnelei  I), 
ci  is  the  a/nuutli  of  the  seatiering  plane,  and 
<’  jh  the  polai  sc  altering  angle. 


The  computer-generated  Mie  intensity  functions  ii(x,^,tj)  and  i2(x,0,tj)  correspond, 
respectively,  to  light  scattered  with  its  electric  field  oriented  perpendicular  and  parallel  to  the 
scattering  plane  (the  "r,  z  plane). 

In  applications  of  Mie  scattering  one  usually  must  consider  scattering  by  a  large  number 
of  particles  (within  a  specified,  finite  volume)  which  may  have  a  multiplicity  of  sizes 
(polydispersion)  and  which  are  viewed  by  a  detector  whose  optics  subtend  a  finite  solid 
angle.  In  such  cases,  Eq.  (1)  must  be  replaced  by 
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where  <f>0  is  the  incident  power  of  the  illuminating  beam;  dn  (f,x)/dx  is  the  number  density 
per  unit  size  interval  of  particles  whose  size  parameter  is  x  and  which  are  located  at  position 
f;  and  4  is  the  projection  of  the  scattering  volume  along  the  optical  axis.  The  differential 
scattering  cross  sections  for  the  two  orthogonal  polarization  states  of  the  scattered  light  are 
related  to  the  Mie  functions  by 

=  —  {x,  6 ,  17)  sin  2  0 

(3) 


=  - (x,  6,  rj)  cos2  <f> 

k2 


and 


da  . 


dn 


(x,  0,  4> .  rf! 


da2 

—  (x,  6,  6,  r,) 


The  subtended  angles  appearing  in  the  limits  of  the  angular  integrals  in  Eq.  (2)  are  obtained 
from  the  standard  geometrical  optics  equations. 

2-3  MIE  SCATTERING  FROM  EXTENDED  REGIONS 

Once  the  scattering  equations  for  localized  particles  (including  polydisperse  samples) 
have  been  obtained,  the  next  task  in  the  feasibility  calculation  is  to  include  effects  caused  by 
the  extension  of  the  source,  as  shown  in  Fig.  1,  for  the  application  at  hand.  The  axial  region 
of  interest  was  approximately  30  cm  (12  in.),  extending  downstream  from  the  source  exit 
plane.  As  shown  in  Fig.  1,  a  thin  slab  of  the  particulate  flow  field  is  irradiated  by  the  laser 
beam.  The  x  axis  is  along  the  flow,  and  the  intensifier  is  located  in  the  y,  z  plane  at  a  distance 
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L  from  the  flow  centerline  and  an  angle  a  from  the  direction  of  the  beam  propagation.  The 
particle  plume  is  illuminated  from  below  by  a  beam-expanded  laser  situated  at  -zL  whose 
beam  has  been  expanded  to  a  divergence  angle  ft,  =  tan'1  Xo/zL  where  xQ  is  half  the  axial 
“look  distance.”  For  these  calculations  x0  was  assumed  to  be  15  cm  (s  6  in.)  and  the  laser 
distance  zL  was  taken  to  be  67.5  cm  (=  27  in.).  Figure  1  shows  scattering  from  a  particle 
which  is  located  at  xp,  yp,  zp,  with  -d/2  <  yp  <  d/2  (where  d  is  the  small  dimension  of  the 
laser  beam)  and  which  is  illuminated  by  an  incident  ray  offset  from  the  laser  axis  by  angle  j3, 
with  -ft,  <  /3  <  ft,.  Note  that  the  plane  of  polarization  of  the  incident  light  is  the  x,  z  plane 
with  respect  to  the  particle  shown,  the  scattering  coordinates  are  r ' ,  9' ,  <j>' ,  and  the  results 
must,  of  course,  be  integrated  over  all  particle  locations  which  can  be  seen  by  the  intensifier. 

The  intensifier  is  composed  of  many  small  resolution  elements  (pixels),  each  of  which 
subtends  a  different  solid  angle  and  looks  at  a  different  region  of  the  flow.  Because  the 
subtended  angle  of  each  pixel  is  very  small,  we  assume  the  light  flux  into  each  pixel  is  the 
same  as  that  at  the  center  of  the  intensifier.  For  simplicity  we  also  ignore  the  spatial 
nonuniformities  of  the  particle  plume  and  of  the  expanded  incident  laser  beam,  assume 
spherical  particles  of  knowm  refractive  index,  and  consider  a  single  particle  size,  or  a 
monodisperse  distribution  function.  Consequences  of  departure  from  these  and  other 
idealizations  will  be  discussed  in  a  later  section. 

Invoking  these  simplifying  assumptions,  we  write  the  total  scattered  flux  received  by  each 
pixel  of  the  intensifier  as 


(4) 


where 


I  =  —  8,  tj'J  sm“  0  +  i  j  0,  cos2  <j>  J  dr  sin  0  "  <\0  '  do ' 


with  the  appropriate  limits  on  r ' ,  8 ’ ,  d> ' ;  A0  is  the  area  of  a  pixel.  The  integration  limits  were 
obtained  as  follows: 

a)  From  Fig.  1  one  sees  that  L  sin  a  -  (d/2)  <  r'  sin  6'  sin  <  L  sin  a  +  (d/2). 

b)  Since  the  intensifier  diameter  is  2.3  cm  (0.9  in.)  and  its  assumed  “look  distance” 
is  30  cm,  its  demagnification  factor  is  M  =  2.3/30  =  0.077. 
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c)  On  the  assumption  of  a  lens  of  focal  length  f  =  6  cm,  the  appropriate  distance 
from  the  intensifier  to  the  plume  centerline  required  to  fill  the  intensifier  face 
was  calculated  with  the  relation 


L  1 

-  =  1  +  - 
f  M 


and  L  was  found  to  equal  86  cm. 
d)  Assuming  symmetry  about  the  y  axis,  we  obtain 


and 


^ —  sin  <*  cos  /3 5  tan  (f> '  1  ^ —  sin  «  cos  £ ^ 

2 

(Lsln“)  +(x0sec/3o) 


<  6  <  tan 


1  f 


’K 


% 


(l  cos  «  +  sec  see  /3^  sec  |l  4 - ^ 

V  7’lJ 


where  Rb,  the  radius  of  the  plume  boundary  at  the  intensifier  axial  station,  was 
assumed  to  be  11.25  cm. 

Using  these  relations  we  determined  the  following  angular  limits: 

a)  75.4°  1  <f>'  <  104.6°  and  45°  <  6'  <  51.2° 
for  oc  =  45  deg 


and 


b)  79.7°  5  <j>'  <  100.3°  and  90°  <  6'  <  97. 5° 

for  oc  =  90  deg 

These  values  were  used  to  calculate  the  integral  I  for  a  range  of  spherical  particle  diameters 
in  the  range  of  I  Atm,  which  for  X  =  514.5  nm  yields  a  size  parameter  x  =  6.11.  In  addition, 
the  index  of  refraction,  17,  was  assumed  to  be  2.68  (Ref.  4).  The  results  are  shown  in  Table  1 . 
Because  of  the  integration  of  the  Mie  functions  over  a  wide  angular  range,  which  means  that 
many  Mie  resonances  are  included  in  the  integral,  no  simple  trend  is  observed  in  Table  1. 
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However,  it  is  seen  that  an  order-of-magnitude  estimate  is  possible  for  each  value  of  a  (e.g., 
IF  as  1  to  10  for  a  =  90  deg). 

Further,  for  the  conditions  of  this  experiment,  the  quantity 

( — - —  — \  =  2.5  x  10'18  citi3 

pe>m  ky 

and  a  typical  number  density  is  assumed  to  be  106  cm  3;  hence, 

$/  =  2.5  x  (l0‘12  to  ltr11) 

/<!>  v  ' 

Finally,  assuming  the  injected  laser  power  to  be  on  the  order  of  one  watt,  it  is  found  that,  for 
these  conditions,  each  intensifier  pixel  receives  a  flux  of  order  107  photons/sec. 

To  determine  the  detectability  of  such  a  photon  rate  per  pixel,  one  must  determine  the 
equivalent  photon  rate  for  the  dark  current  of  the  tube.  Assuming  the  dark  current  density 
to  be  equivalent  to  an  input  photon  flux  density  of  5  x  104  photons/ sec/ cm2*  and  taking  the 
resolution  element  to  be  0.07  mm,  one  finds  the  photon  noise  rate  of  a  single  resolution 
element  to  be  approximately  two  photons/sec;  obviously,  the  single-to-noise  ratio  ( *=  106  to 
107)  is  adequate  for  this  study. 

An  additional  intensifier  system  parameter  of  interest  is  the  photon  gain  (g),  i.e.,  the 
number  of  photons  at  the  output  phosphor  per  photon  at  the  photocathode.  The  image 
intensifier  parameters  which  are  appropriate  for  the  device  used  for  this  study  are 

Electron  Gain  (G)  =  500  output  electrons/photocathode  electron 

Photon  Conversion  Efficiency  (E)  =  0.05  output  photons/output  electron/eV 

Effective  Accelerating  Voltage  (V)  =  2.5  keV 

Photocathode  Efficiency  (f)  at  500  nm  =  0.10 

The  photon  gain  is  easily  found  to  be 

g  . V . L ■ G 

or 

=  6  \  10  oulpui  photon*-  inpuL  photon. 


*1TT  F-4112  intensifier  tube 
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Therefore,  assuming  an  optical  system  transmission  factor  of  approximately  0.15,  one  finds 
the  effective  photon  gain  of  the  intensifier  to  be  on  the  order  of  103. 

To  compare  such  a  system  to  straight  photographic  recording,  it  is  easiest  to  suppose  the 
intensifier  output  to  be  recorded  by  a  high-speed  cine  camera  whose  framing  rate  defines  the 
time  resolution  of  the  measurement.  If  such  a  system  is  compared  with  detection  by  this 
same  cine  camera  alone,  it  is  obvious  that  the  intensifier  system  possesses  an  inherent 
advantage  given  by  the  photon  gain,  which,  in  this  case,  was  on  the  order  of  103.  However, 
one  can  obtain  a  more  decisive  comparison  by  considering  the  minimum  detectable  signal 
for  photographic  recording.  A  typical  minimum  energy  density  for  film  exposure  is  on  the 
order  of  0. 1  erg/cm2.  For  an  area  corresponding  to  that  of  the  minimum  resolution  element 
of  the  intensifier  one  finds  that  approximately  106  photons  at  500  nm  are  required  to  exceed 
the  typical  film  inertia.  For  the  particulate-laden  flow  discussed  in  a  previous  paragraph,  it 
was  found  that  the  anticipated  Mie  scattered  signal  per  pixel  was  on  the  order  of  107 
photons/sec.  Therefore,  for  a  0.  l-sec  duration  flow  field,  approximately  106  photons  are 
scattered  onto  a  pixel  element  area.  Consequently,  the  detectability  of  such  Mie  scattering  by 
photographic  methods  alone  is  quite  marginal  for  exposure  times  significantly  less  than  0.1 
sec. 


3.0  EXPERIMENTAL  APPARATUS 

On  the  basis  of  the  results  of  the  feasibility  calculation  of  the  previous  section  the  image 
intensifier  system,  sketched  in  Fig.  1,  was  used  to  visualize  the  particle-laden  flow  field, 
which  was  immersed  in  a  nominal  Mach  13  flow.  The  particle  injection  source  was  a  blunt 
cone  with  a  base-mounted  nozzle  of  approximately  1-mm  diameter.  The  nozzle  centerline 
was  directed  outward  at  45  deg  from  the  model  axis,  and  the  particle  exhaust  was  comprised 
of  nominal  l-/i-diam  particles  borne  by  helium  (He)  carrier  gas.  Further  details  of  the  flow 
fseld  and  particle  source,  although  of  importance  for  the  analysis  of  the  particle/ flow-field 
interaction,  will  not  be  discussed  in  detail  since  the  purpose  of  this  report  is  primarily  to 
describe  the  visualization  technique  itself.  Flow-field  properties  will  be  included  only  when 
they  are  essential  for  evaluating  this  diagnostic  technique. 

The  experimental  geometry  utilized  in  the  visualization  measurements  and  sketched  in 
Fig.  1  is  shown  in  more  detail  in  Fig.  3.  Source  illumination  for  the  particulate  plume  was 
supplied  by  a  Coherent  Radiation  Model  CR-3  argon-ion  laser  operating  in  the  continuous 
mode  at  a  wavelength  of  514.5  nm  with  a  nominal  power  level  of  1.5  watts.  The  laser  (LA) 
was  mounted  beneath  the  tunnel  on  a  vibration  isolation  table,  and  it  produced  a  coherent, 
monochromatic  beam  of  light  approximately  2  mm  in  diameter.  This  beam  was  passed 
through  the  blades  of  a  chopper  (c),  which  modulated  the  laser  .beam  at  100  Hz,  and  was 
then  directed  into  an  optical  system  by  a  dielectric  mirror. 
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The  optical  system  consisted  of  a  spatial  filter  (SF),  cylindrical  lens  (LN),  and  two 
limiting  apertures  which  produced  a  fan-shaped,  extended  beam  approximately  2  mm  thick 
which  was  directed  through  the  225-deg  port  into  the  wake  flow  region  of  the  model  and 
passed  through  a  plane  formed  by  the  tunnel/model  centerline  and  the  axis  of  the  particulate 
ejection  nozzle. 

The  beam  was  essentially  uniform  in  intensity  across  its  entire  spatial  extent  and  was  19.1 
cm  (7.55  in.)  in  width  at  the  tunnel  centerline  and  20.3  cm  (8.0  in.)  in  width  at  the  particulate 
ejection  orifice  (see  Fig.  4).  The  expanding  laser  beam  exited  the  test  section  through  a 
quartz  window  and  was  intersected  by  six  transmission  monitors  located  at  the  45-deg  port. 

The  image  intensifier  (lT)/fast  camera  (FC)  detection  system,  shown  in  Fig.  4,  was 
positioned  so  that  the  viewing  angle,  a,  was  90  deg  with  respect  to  the  plane  of  incident  laser 
light.  Since  only  the  particles  passing  through  the  laser-illuminated  region  contributed  to  the 
scattered  light  (for  the  case  of  visible  scattered  radiation)  the  photographic  images  of  the 
plume  viewed  in  this  manner  represented  a  two-dimensional  slice  through  the  plume. 

The  detection  system  consisted  of  a  collection  lens  (LI)  which  gathered  the  scattered 
radiation  and  focused  the  image  of  the  visualization  area  onto  the  active  surface  of  the  image 
intensifier.  The  intensifier  used  was  an  ITT  F-41 12  with  an  S-20R  spectral  response,  which  is 
shown  in  Fig.  5.  The  electron  gain  of  the  tube  is  nominally  500  for  an  operating  voltage  of 
5,000  V.  The  output  of  the  intensifier  was  focused  by  a  second  lens  (L2)  onto  the  film  of  a 
high-speed,  16-mm  Mitchell  Monitor  Camera.  Filters  F  provide  spectral  isolation  of  the 
scattered  radiation.  The  use  of  several  different  filter  combinations  provided  a  variety  of 
visualization  options.  With  the  incident  laser  beam  unblocked  by  the  chopper  and  a 
514.5-nm  narrow  bandpass  filter,  only  Mie  scattered  radiation  from  the  particles  was 
recorded.  This  filter  had  a  bandwidth  of  1  nm  full  width  at  half  intensity  maximum. 
Utilizing  a  near-infrared  bandpass  filter,  which  transmitted  in  the  region  from  700  nm  to  1.5 
/tm,  in  conjunction  with  the  laser  beam  blocked  by  the  chopper  blades  provided  for  an 
investigation  of  the  near-infrared  particle  luminosity.  A  combination  of  the  data  obtained 
by  using  both  filters  provided  For  transmission  in  the  regions  514.0  to  515.0  nm  and  700  nm 
to  1.5  and  produced  a  signal  representative  of  the  combination  of  near-infrared 
particle  luminosity  and  Mie  scattered  radiation.  With  the  chopper  cycling  the  laser  source 
radiation  from  on  to  off,  all  of  the  above  visualization  cases  were  obtained. 

The  intensifier  output  was  recorded  using  a  Mitchell  camera  which  was  operated  at  a  rate 
of  400  frames  per  second,  resulting  in  a  time  between  frames  of  2.5  msec;  the  frame  or 
photographic  time  resolution  is  less  than  2.5  msec  and  is  estimated  to  be  approximately  1 
msec.  Reference  times  were  indicated  by  comparison  with  the  tunnel  arc  voltage.  Selected 
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film  data  frames  were  digitized  utilizing  the  VKF  image/ Analyzer  Digitizer  System  (IADS) 
and  stored  on  computer  disk  file  for  subsequent  data  processing.  (The  IADS  is  a  PDP-11- 
controlled  film  densitometer  device;  it  measured  and  stored  the  data  film  density  as  a 
function  of  spatial  position.)  Subsequent  data  reduction  included  calibration  data  for  the 
effects  of  laser  beam  intensity  profile  variations,  image  intensifier  response,  and  the  addition 
of  spatial  reference  grids.  Computer-based  data  reduction  routines  were  used  in  analyzing 
data  film  frames  with  respect  to  light  intensity  level,  and  final  photographic  records  were 
produced  with  grayness  (or  color)  levels  representing  scattered  intensity  levels  and  lines  of 
constant  light  intensity.  Computer  routines  also  provided  for  intensity  subtraction  on  a 
spatial  point-by-point  basis  so  that  successive  frames  of  near-infrared  combined  with  Mie 
scattered  radiation  frames  could  be  combined  with  alternating  pure  infrared  frames  to 
produce  the  isolated  Mie  data. 

Six  calibrated  silicon  photodetectors  were  used  to  measure  the  degree  of  attenuation  of 
the  laser  radiation  as  it  passed  through  the  particulate  plume. 

These  measurements  were  made  at  several  spatial  locations  within  the  particle  plume  as 
indicated  in  Fig.  6.  The  laser  transmission  path  to  each  detector  subtended  a  small  volume  of 
the  plume,  and  the  attenuation  of  the  laser  power  along  each  individual  path  was  recorded  as 
a  function  of  time  by  the  Tunnel  F  Transient  Data  System  (TDS). 

Just  prior  to  each  tunnel  firing,  a  laser  transmission  calibration  at  tunnel  vacuum  was 
recorded  on  the  TDS  to  provide  a  maximum  transmission  (no  particle)  reference.  During  the 
actual  run,  the  TDS  recorded  the  laser  transmission  channels  at  1-msec  intervals. 
Transmission  factors  were  then  determined  at  selected  times  for  further  data  analysis.  Along 
the  laser  beam  optical  path,  s,  energy  is  lost  from  the  incident  beam  by  scattering  and 
absorption.  In  the  absence  of  multiple  scattering,  the  energy  flux  of  transmitted  light  is 


^TR  =  °oexP 


(5) 


where  uexl(x,s)  is  the  size  parameter  (x)-  and  position  (s)-dependent  extinction  cross  section 
and  S  is  the  total  path  length.  Assuming  the  distribution  function  to  be 

<?n  (x,  s) 
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where  n  is  the  uniform  particle  density,  if  the  distribution  is  both  monodisperse  and  spatially 
uniform,  one  finds  that 

®'l  u  =  4>o  exP  [- "  ‘'cm  (*o)  S]  ^ 


Therefore,  if  the  size  parameter,  x0,  for  such  a  sample  is  known,  aext  is  calculable. 
Furthermore,  the  path  length  (S)  can  be  determined  from  the  visualization  photographs,  and 
the  transmission  ratio  r  =  can  be  measured.  With  these  parameters  in  hand  the 

particle  density  along  the  laser  beam  path  is  determined. 

Additionally,  the  mass  flow  rate  (mp)  of  the  particles  can  be  related  to  the  attenuation  of 
the  laser  radiation  by  the  particle  plume  in  the  following  manner: 


where  is  the  density  of  the  solid  particulate,  vp  is  the  particulate  flow  speed  (flow  speed),  f 
is  the  length  of  the  intersected  particulate  plume,  D  is  the  mean  particle  diameter,  trD  is  the 
mean  extinction  cross  section,  and  r  is  the  transmission  fraction  determined  as  described 
above.  The  known  optical  properties  of  the  particles  and  the  mean  particle  diameters  of  the 
initial  samples,  determined  through  independent  analysis,  were  used  in  conjunction  with 
measured  transmission  fractions  to  compute  experimentally  determined  particulate  mass 
flow  rates.  Transmission  path  lengths  for  these  computations  were  determined  from  the 
visible  visualization  photographs.  Results  could  then  be  compared  with  mass  flow  rates 
determined  from  model  canister  weight  measurements  and  mass  flow  times. 

4.0  RESULTS 


4.1  VISUALIZATION  RESULTS 

Exemplary  flow  visualization  results  are  presented  in  Figs.  7  through  13  and  summarized 
in  Table  2.  Data  for  several  different  run  conditions  and  particle  compositions  are  presented 
to  illustrate  the  wide  range  of  possible  visualization  data  obtainable  under  the  techniques 
used  in  this  study. 

Data  from  Fig.  7  consisted  of  Mie  scattered  radiation  combined  with  near-infrared 
particulate  luminosity.  In  Fig.  8  pure  infrared  signal  was  recorded  with  the  laser  source 
radiation  blocked  by  the  chopper  blade.  Here  the  particulate  plume  emanates  from  the 
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nozzle,  located  on  the  left  side  of  Fig.  8,  and  expands  to  the  right.  Flow  direction  is  from  left 
to  right  for  all  data  figures.  As  can  be  clearly  seen,  the  infrared  signal  begins  somewhat 
downstream  from  the  nozzle  exit  and  in  a  region  paralleling  the  bow  shock  from  the  model. 
Indications  are  that  the  particles  require  a  short  distance  or  time  to  interact  with  the  flow 
before  the  particle  temperature  rises  sufficiently  to  produce  an  infrared  signal.  Figure  7  was 
originally  recorded  with  the  laser  beam  on  as  a  combination  of  visible  and  near-infrared 
radiation.  In  order  to  obtain  a  pure  visible  signal,  the  pure  IR  signal  V  (intensity)  from  Fig.  8 
was  subtracted  point  by  point  from  Fig.  7.  One  can  see  from  Fig.  7,  then,  that  the  visible 
scattered  radiation  from  the  particles  clearly  traces  the  particles’  exit  From  the  nozzle,  their 
interaction  with  the  bow  shock,  and  their  propagation  downstream.  Obviously,  absolute 
signal  levels  have  been  modified  in  these  computer  representations  for  the  sake  of  clarity. 
The  overlying  grid  lines  represent  1-in.  squares  within  the  flow  and  are  provided  for  spatial 
reference.  The  particles  were  pure  tungsten  (W). 

Figure  9  presents  typical  data  from  a  tungsten  carbide  particulate  plume.  Here  the  data 
represent  almost  entirely  visible  scattered  radiation.  Again,  the  particles  leave  (he  nozzle  in 
high  concentration  and  expand  as  they  move  downstream  until  they  are  dispersed  to  the 
point  that  signal  levels  are  below  threshold. 

Figure  10  illustrates  a  case  of  scattering  from  silicon  carbide  (SiC)  particles.  The 
particulate  temperature  for  this  run  was  very  low  in  comparison  with  prior  runs,  and  an 
absence  of  near-infrared  radiation  can  be  discerned. 

Figures  II  and  12  graphically  illustrate  particulate  heating  upon  interaction  with  the 
tunnel  flow.  In  Fig.  1 1 ,  very  early  in  the  run,  the  tungsten  carbide  (WC)  particles  display  a 
very  weak  IR  signal  and  a  highly  dispersed  visible  signal.  However,  much  later  in  the  run 
(Fig.  12),  particulate  temperature  has  increased,  and  a  very  strong  infrared  signal  is 
observed.  Considerably  different  run  conditions  from  previous  data  frames  were  present,  as 
can  be  seen  from  the  photographic  data  summary. 

Finally,  Fig.  13  illustrates  a  very  hot  particulate  (1.578K)  plume,  and  a  very  strong 
infrared  signal  is  observed. 

4.2  TRANSMISSION  RESULTS 

During  the  course  of  this  measurement  program  the  average  particulate  mass  flow  rate 
(mp)  was  experimentally  determined  by  using  the  mass  loss  of  the  particulate  generator  and 
the  total  ejection  time  of  the  particles.  It  is  recalled  from  Eqs.  (6)  and  (7)  that  if  the  particle 
distribution  is  uniform  in  space  and  monodisperse  in  size,  measurement  of  the  transmission 
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factor  (t)  yields  the  particle  mass  flow  rate.  Even  though  the  Mie  scattering  photographs 
vividly  show  a  violation  of  the  assumption  of  spatial  uniformity,  it  is  of  interest  to  apply  Eq. 
(7)  and  the  measured  values  of  r  to  determine  rhp.  The  results  of  this  exercise  are  shown  in 
Table  3  for  four  different  tunnel  runs.  From  Table  2  one  can  see  that  order-of-magnitude 
agreement  is  obtained  in  all  cases  with  the  mass-loss  measurement  of  mp. 

5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  this  report  was  to  demonstrate  the  capabilities  of  a  Mie  scattering-based 
image  intensifier  system  for  particulate  diagnostics,  and  it  is  suggested  that  this  goal  has 
been  achieved.  The  sensitivity  of  this  system,  which  employed  a  high-speed  cine  camera  for 
recording,  obviously  exceeds  that  of  the  high-speed  camera  alone,  and  the  capability  for 
qualitative  diagnostics  of  a  broad  spatial  field  has  been  demonstrated.  Although  it  is  perhaps 
obvious,  it  is  significant  to  note  that  this  diagnostics  capability  does  provide  spatial 
resolution;  i.e.,  scattering  occurs  only  in  the  region  of  iaser  illumination  which,  in  this  case, 
approximates  a  plane.  Further,  the  time  resolution  of  the  system,  which  was  on  the  order  of 
1  msec,  was  limited  by  the  recording  camera  speed  rather  than  by  the  sensitivity  of  the 
intensifier,  which,  it  should  be  noted,  is  a  comparatively  low  gain  tube. 

Quantitative  results  of  the  particulate  spatial  profiles  are  obtainable  only  if  certain 
assumptions  are  valid.  Included  in  these  assumptions  are  knowledge  of  the  size  distribution 
function,  including  a  monodisperse  distribution,  and  the  index  of  refraction  of  the  particles. 
Further,  use  of  Eq.  (7)  to  relate  the  transmission  function  to  the  particulate  mass  flow 
requires  an  even  more  stringent  assumption  of  spatial  uniformity  of  the  particles,  which 
Figs.  7  through  13  show  would  obviously  not  be  valid.  Consequently,  agreement  of 
particulate  mass  flux  rates  as  determined  by  mass  loss  and  transmission  techniques  to  within 
no  better  than  an  order  of  magnitude  is  not  surprising;  in  fact,  when  one  acknowledges  not 
only  the  spatial  nonuniformity  of  the  particles  but  also  their  known  polydispersity,  the 
order-of-magnitude  agreement  shown  in  Table  3  is  quite  gratifying.  Quite  clearly,  for  a 
spatially  nonuniform  flow  field  a  local  measurement  of  particle  speed  and  number  density  is 
required  for  a  local  mass  flow-rate  determination.  Such  measurements  of  flow  speed  have 
been  performed  and  will  be  the  subject  of  a  future  report.  Further,  for  a  monodisperse 
distribution  of  injected  particles  the  local  number  density  can  be  determined  in  a  straight¬ 
forward  manner  using  an  in  situ  calibration.  For  a  polydisperse  distribution  either  particle 
sizing  must  be  performed  or  it  must  be  assumed  that  the  size  distribution  function  is  spatially 
invariant;  in  either  event,  the  local  particulate  mass  flow  rate  can  be  determined.  In 
addition,  the  isointensity  contours  of  Figs.  7  through  13  can  be  reliably  and  quantitatively 
interpreted  as  particulate  isodensity  contours. 
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In  conclusion,  when  one  incorporates  the  just-mentioned  calibration  procedure  into  the 
intensifier  detection  and  analysis  system,  the  technique  represents  a  significant  enhancement 
in  diagnostic  capabilities  for  study  of  particle-flow  field  interactions. 
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re  1.  Particle-laden  flow-field  configuration. 
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Figure  2.  Scattering  geometry  for  single  particle  with  light  incident 
along  Z  axis  and  polarization  vector  along  X  axis. 
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View  Looking  Upstream 


Nomenclature 


Chopper 

Optical  Filter 

Fast  Camera  (Camera  No*  1) 

Image  In tens i tier 

Laser 

Lens 

Spatial  Filter 
Transmission  Monitor 
Viewing  Angle 


Figure  3.  Experimental  configuration:  Flow  visualization  and 
transmission  measurements. 
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Relative  Sensitivity,  raa 


Figure  5.  Typical  intensifier  response  curve. 
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A  EDC-TR  -79-42 


PD  -  Photodetector 
PD6  PD5  PD4  PD3  PD2  PD1 


Laser 


Figure  6.  Transmission  channel  configuration. 
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Figure  7.  Flow  visualization  of  tungsten  particulate  flow.  Run  5,986, 
visualization  signal  source  resulting  from  Mie  scattered 
radiation  combined  with  near-infrared  particulate  luminosity. 
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to 


Figure  8.  Flow  visualization  of  tungsten  particulate  flow.  Run  5,986, 
visualization  signal  source  consisting  entirely  of  near-infrared 
particulate  luminosity. 
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Figure  9.  Flow  visualization  of  tungsten  carbide  particulate  flow, 
Run  5,990,  visualization  signal  source  consisting  entirely 
of  Mie-scattered  radiation. 


A  E  DC-TR -79-42 


Figure  10.  Flow  visualization  of  silicon  carbide  particulate  flow. 
Run  5,991,  visualization  signal  source  resulting  from 
Mie-scattered  radiation  combined  with  near-infrared 
particulate  luminosity. 
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Run  5,996,  Frame  46 


Figure  11.  Flow  visualization  of  tungsten  carbide  particulate  flow. 

Run  5,996  at  1 1 1  msec  from  tunnel  arc  voltage,  visualization 
signal  source  consisting  of  very  weak  near-infrared  particle 
luminosity  and  Mie-scattered  radiation. 
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Figure  12.  Flow  visualization  of  tungsten  carbide  particulate  flow. 

Run  5,996  at  131  msec  from  arc  voltage,  visualization  signal 
source  consisting  of  strong  near-infrared  particulate 
luminosity  and  Mie-scattered  radiation. 


AEDC-TR-79-42 


Figure  13.  Flow  visualization  of  tungsten  particulate  flow.  Run  5,997, 

with  a  particle  plume  temperature  of  approximately  1,578  K, 
visualization  signal  source  consisting  of  near-infrared 
particulate  luminosity  and  Mie-scattered  radiation. 
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Table  1.  Flux  Integral  IF 


X 

a  ®  45  def* 

a  a  90  deft 

5,11 

31.9 

7.9 

6.11 

57.6 

2.3 

7.11 

146.9 

3.8 

8.11 

38.4 

26.4 

Table  2.  Exemplary  Data  Photograph  Summary 


Number ' 

Particle 

Time 
from  Arc 
Voltage,  msec 

Re  x 

4) 

10"6/ft 

T0,  K 

Equivalent 
Blackbody 
Plume  Temp . ,  K 

Radiation  Source 

7 

V 

96 

2.88 

2,042 

1,302 

Visible  (514.5  nm) 

Near-IR  (700  nm-1,5  um) 

8 

w 

101 

**— 

— 

— 

Near-IR  (700  nm-1.5  \im) 

9 

wc 

89 

2.25 

2,210 

1,138 

Visible  +  (Weak)  IR 

10 

SiC 

111 

2.74 

1,964 

-Lov- 

Visible  +  (Weak)  IR 

11 

WC 

111 

0.8 

2,232 

1,113 

Visible  +  (Very  Weak)  IR 

12 

WC 

131 

— 

— 

— 

(Strong)  IR  +  Visible 

13  . 

w 

84 

0.47 

2,868 

1,578 

(Very  Strong)  IR  +  Visible 

fable  3.  Mass  Flow  Rate  —  Exemplary  Data  Summary 


Run 

Particle 

Ap  (Transmission) , 
ftm/sec 

(Weight), 

Ert/aec 

1 

SiC 

4.22 

1.0 

2 

W 

3,50 

4.4 

3 

SiC 

1.91 

1.33 

4 

W 

3.91 

4.3 

33 


Table  4.  Exemplary  Transmission  Data 


Figure 

Number 

Time  from 
Arc  Voltage, 
msec 

Chopper 

Particle 

Material 

PD1 

Particle  Concentration  (Particles/cm^) 
PD2  PD3  PD4 

PD5 

PD6 

7 

124 

Ye3 

W 

— 

- 

2.63 

X 

106 

— 

_ .. 

_ 

8 

116 

Yes 

wc 

1.51 

X 

105 

1.48 

X 

105 

1.31  x  1<>5 

6.88  x 

104 

6.11 

X 

104 

5.27 

x  104 

9 

ill 

Yes 

SiC 

2.76 

X 

105 

2.37 

X 

105 

— 

1 . 18  x 

105 

5.98 

X 

104 

2.15 

x  104 

10 

111 

No 

W 

5*07 

X 

105 

8.16 

X 

105 

A. 31  x  105 

2.44  x 

105 

3.36 

X 

105 

3.04 

x  105 

il 

120 

No 

SiC 

2.97 

X 

105 

2.69 

X 

105 

1.07  x  105 

4.28  x 

104 

3.11 

X 

104 

5*97 

x  104 

12 

125 

No 

wc 

6.22 

X 

104 

9.89 

X 

104 

7.33  x  104 

5.91  x 

104 

5.27 

X 

104 

2.43 

x  104 

13 

119 

No 

w 

1.44 

X 

105 

2.57 

X 

105 

2.86  x  105 

1.50  x 

105 

1.03 

X 

105 

6.91 

x  104 

PD  -  Photodetector 


AEDC-TR- 79-42 


Abeam 

Ao 

D 

d 

E 

F 

FC 

f,  LI,  L2 
G 

g 

I 

lo 

IF 

II 

i|,2 

k 

L 

LN 

i 

M 

nip 

n 

n 


NOMENCLATURE 

Area  of  source  beam 

Area  of  pixel 

Mean  particle  diameter 

Thickness  of  incident  laser  sheet 

Photon  conversion  efficiency 

Optical  filters  within  the  collection  optics  path 

Fast  camera 

Focal  lengths  of  lenses  within  the  image  intensifier  optical  path 
Electron  gain 

Image  intensifier  photon  gain 

Intensity  of  scattered  radiation 

Incident  laser  intensity 

Flux  integral 

Image  intensifier 

Mie  intensity  functions 

Wave  number  of  incident  light 

Distance  from  flow  centerline  to  image  intensifier 

Cylindrical  lens 

Position  of  scattering  center  along  laser  path  (Section  2);  length  of 
intersected  particulate  plume  (Section  3) 

Demagnification  factor 

Mass  flow  rate  of  particles 

Particle  number  density 

Particle  number  density  averaged  over  path  length,  s 
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Rb  Radius  of  plume  boundary 

r  Distance  from  scattering  event  to  detector 

r,#,0  Polar  coordinates  (axis  shown  in  Fig.  2) 

5  Total  path  length  of  laser  beam 

SF  Spatial  filter 

s  Laser  beam  optical  path  along  which  scattering  and  absorption  occur 

V  Effective  accelerating  voltage 

vp  Particulate  flow  velocity 

X  Particle  size  parameter 

x,y,z  Rectangular  coordinate  system  as  shown  in  Fig.  1 

x0  Half  the  axial  “look  distance" 

Zl  Distance  from  flow  centerline  to  laser  source 

a  Angle  separating  directions  of  laser  propagation  and  intensifier  optical  axis 

(Section  2) 

0  Angle  of  divergence  of  injection  beam 

6  Beam  divergence  angle  with  -/30  <  /?  <  /?0 

f  Photocathode  efficiency 

rj  Index  of  refraction  relative  to  the  medium 

Vi  Real  part  of  index  of  refraction 

0  Polar  scattering  angle  (Section  2) 

k  Imaginary  part  of  index  of  refraction 

X  Wavelength  of  incident  light  in  the  scattering  medium 

P s  Density  of  solid  particulate 

o  Scattering  cross  section  integrated  over  solid  angle  subtended  by  detector 

^xt  Extinction  cross  section 
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<*0 

T 

$ 

SUBSCRIPTS 

o 

1 

2 

s 

TR 


Mean  extinction  cross  section 
Transmission  ratio 
Scattered  power 
Incident  laser  power 
Azimuth  of  the  scattering  plane 


incident 

Denotes  scattering  state  with  electric  vector  perpendicular  to  scattering 
plume 

Denotes  scattering  state  with  electric  vector  parallel  to  scattering  plane 
Denotes  parameters  specific  to  scattering  volume  and  angles 
Denotes  transmitted  beam 
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